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Abstract 



The present study deals with spatially homogeneous and totally anisotropic locally rotationally symmetric 
■ (LRS) Bianchi type I cosmological model with variable G and A in presence of imperfect fluid. To get the 

deterministic model of Universe, we assume that the expansion (8) in the model is proportional to shear (<r). 
This condition leads to A = £B n , where A, B are metric potential. The cosmological constant A is found 
to be decreasing function of time and it approaches a small positive value at late time which is supported 
by recent Supernovae la (SN la) observations. Also it is evident that the distance modulus curve of derived 
I , model matches with observations perfectly. 
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1 Introduction 



Recent astronomical observations of type la supernovae with redshift parameter z < 1 (Perlmutter et al. 1997, 
1998, 1999; Riess et al. 1998, 2004; Garnavich et al. 1998a, 1998b), Wilkison Microwave Anisotropy Probe 
(WMAP) (Spergel et al. 2003) etc. provide evidence that we may live in low mass density Universe i. e. Q ~ 0.3 
(Riess 1986). The predictions of observations lead to a convincing belief in modern cosmology that a part of 
Universe is filled up with dark energy (f2 ~ 0.7), which may be addressed by suitable cosmological constant. 
There are significant observational evidence that the expansion of the Universe is undergoing a late time accel- 
eration (Perlmutter et al. 1997, 1998, 1999; Riess et al. 1998, 2004; Efstathiou et al. 2002; Spergel et al. 2003; 
\Q ; Allen et al. 2004; Sahni and Starobinsky 2000; Peebles and Ratra 2003; Padmanabhan 2003; Lima 2004). This, 
in other words, amounts to saying that in the context of Einstein's general theory of relativity some sort of dark 
energy, constant or that varies only slowly with time and space dominates the current composition of cosmos. 
The origin and nature of such an accelerating field poses a completely open question. The main conclusion of 
these observations is that the expansion of the Universe is accelerating. 

Among many possible alternatives, the simplest and most theoretically appealing possibility for dark energy 
is the energy density stored on the vacuum state of all existing fields in the Universe, i.e., p v = g^y, where A is 
the cosmological constant. However, a constant A cannot explain the huge difference between the cosmological 
constant inferred from observation and the vacuum energy density resulting from quantum field theories. In an 
attempt to solve this problem, variable A was introduced such that A was large in the early universe and then 
decayed with evolution (Dolgov 1983). Since the pioneering work of Dirak (1938), who proposed a theory with a 
time varying gravitational coupling constant G(t), a number of cosmological models with variable G and A have 
been recently studied by several authors (Arbab 2003; Sistero 1991; Sattar and Vishwakarma 1997; Pradhan 
and Chakrabarty 2001; Singh et al. 2008). 

To describe the relativistic theory of viscosity, Eckart (1940) made the first attempt, the theories of dissipa- 
tion in Eckart formulation suffers from serious short- coming, viz., causality and stability (Hiskock and Lindblom 
1985; Hiskock 1986) regardless of the choice of equation of state. The problem arises due to first order nature of 
the theory, since it considers only first order deviation from equilibrium. It has been shown that the problems 
of the relativistic imperfect fluid may be resolved by including higher order deviation terms in the transport 
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equation (Hiskock and Salmonson 1991). Isreal and Stewart (1970) and Pavon (1991) developed a fully relativis- 
tic formulation of the theory taking into account second order deviation terms in the theory, which is termed 
as "transient" or "extended" irreversible thermodynamics (EIT). The crucial difference between the standard 
Eckart and the extended Isreal-Stewart transport equations is that the latter is a differential evolution equations, 
while the former is an algebraic relation. In irreversible thermodynamics, the entropy is no longer conserved, but 
grows according to the second law of thermodynamics. Bulk viscosity arises typically in the mixtures either of 
different species or of species but with different energies. The solution of the full causal theory are well behaved 
for all the times. Therefore the best currently available theory for analyzing dissipative processes in the Universe 
is the Full Isreal-Stewart theory (FIS). Several authors (Kremer et al. 2003; Singh and Beesham 2000; Debnath 
et al. 2007; Singh and Kale 2009) and recently Yadav (2010, 2011) have obtained cosmological models with 
dissipative effects. Pradhan et al. (2004); Singh (2009); Singh and Kumar (2009) and Bali and Kumawat (2008) 
have studied matter filled imperfect fluid in different physical context. 

The simplest of anisotropic models are Bianchi type-I homogeneous models whose spatial sections are flat 
but the expansion or contraction rate are direction dependent. For studying the possible effects of anisotropy in 
the early Universe on present day observations many researchers (Huang 1990; Chimento et al. 1997; Lima 1996; 
Lima and Carvalho 1994; Pradhan and Singh 2004; Pradhan and Pandey 2006; Saha 2005, 2006a, 2006b) have 
investigated Bianchi type-I models from different point of view. In this paper, we present the exact solution of 
Einstein's held equations with variable G and A in LRS Bianchi I space-time in presence of imperfect fluid as a 
source of matter. The paper has following structure. In section 2, the metric and field equations are described. 
The section 3 deals with the exact solution of the field equations and physical behaviour of the model. The 
distance modulus curve is described in section 4. At the end we shall summarize the findings. 



2 The Metric and Field Equations 

We consider the LRS Bianchi type I metric of the form 

ds 2 = -dt 2 + A 2 dx 2 + B 2 (dy 2 + dz 2 ) , (1) 
where, A and B are functions of t only. This ensures that the model is spatially homogeneous. 
The energy-momentum tensor T? for bulk viscous fluid is taken as 

T; = (p + p + nK^ + b + n) 5 ;-, (2) 

where p is the isotropic pressure; p is the energy density of matter; LI is the bulk viscous stress; v % = (0, 0, 0, 1) 
is the four velocity vector satisfying the relations 

vy = -1. (3) 

The bulk viscous stress is given by 

n = -£*4 (4) 

where £ is the bulk viscosity coefficient. 

The Einstein's field equations with cosmological constant may be written as 

/,'; - ^g)R - Ag{ = >-(,! :. (5) 
The Einstein's field equations © for the line-element JTJ) lead to the following system of equations 

2^ + || = -8 7 rG(p + n)+A ! (6) 
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Here, and in what follows, sub in-dices 4 in A, B and elsewhere indicates differentiation with respect to t. 
In view of vanishing divergence of Einstein tensor, we get 



8ttG 



, Aa Ba 

P4 + (p + p + ii) 4 + 2^ 



+ 8npG 4 + A 4 = 0. (9) 



A B 

Using Eq. the energy conservation equation © splits into two equations 

P4 + (p +p + n)(^± + 2^f) =o, (io) 

and 

87rpG 4 + A 4 = 0. (11) 
The average scale factor (a) of LRS Bianchi type I model is defined as 

a={AB 2 )i. (12) 

The spatial volume (V) is given by 

V = a 3 = AB 2 . (13) 
We define the mean Hubble parameter (H) for LRS Bianchi I space-time as 

a 4 1 / Aa Ba \ , , 

ff =i=3-U + 2 i)- < 14 > 
The expansion scalar (8), shear scalar (a) and mean anisotropy parameter (Am) are defined as 

e = 3ff = ^ + 2§, (15) 



(16) 



2— 1 x ' 

3 Solutions of the Field Equations 

The system of eqs. ©— ©, (jTDJ) and (ITT1) is employed to obtain the cosmological solution. The system of 
equations is not closed as it has seven unknown (A, B, p, P n, G and A) to be determined from five equations. 
Therefore, two additional constraint relating these parameter are required to obtain explicit solutions of the 
system. 

Firstly, we assume that the expansion (9) in the model is proportional to the shear (a). This condition leads to 

A = tB n , (18) 
where t and n are constant of integration and positive constant respectively. 

Following Luis (1985), Johari and Desikan (1994), Singh and Beesham (1999) and recently Singh and Kale 
(2009), we assume the well accepted power law relation between gravitational constant G and scale factor a as 



where Go and m are positive constants. 
Equations ©, © and (TSJ lead to 



The solution of equation (|2"0"|) is given by 



G = G a m , (19) 



^ + (n+l)§=0. (20) 



B= (fat + ko)^, (21) 
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where fco and fci are the constants of integration. 
From equations (p~8|) and (|2T|) . we obtain 



A = £(kit + k )^+t . 

The rate of expansion in the direction of x, y and z are given by 

A 4 nki 1 
ti x — — r — 



A (n + 2) (fa* + fco) 

= 7f - fcl 1 
a 2 (n + 2)(M + fco)' 

The mean Hubble's parameter (H), expansion scalar (9) and shear scalar (a) are given by 

fci 



H 



3(fcit + fc )' 
fci 

(fcit + fco)' 



a 2 = 



(n- l) 2 fc 2 



3(n + 2) 2 (M + fc ) 2 ' 

The spatial volume (V), mean anisotropy parameter (Am), and average scale factor (a) are found to be 

V = £(fci< + fc ), 
2(n- l) 2 



From equations ([26]) and (|27|) . we obtain 



Am ~ (n + 2)2 ' 

a- _ (n-1) 
~ V3(n + 2) ' 

a= [^(M + fco)]*- 

For specification of £, we assume that the fluid obeys the equation of state of the form 

P = IP, 

where 7(0<7<l)isa constant and it is termed as Equation of state parameter (EoS parameter). 
Differentiating equation ©, we obtain 



8ttGp4 + 8ttG 4 p + A 4 



(2n + l)k\ 



(n + 2) 2 (k 1 t + k ¥ 



From equation dTTJ, 1(15)1. ((51)) and ((22)), we obtain 



3(2n + l)fci 



8nG £T>m(n + 2) 2 (k 1 t + fco)^ 



(22) 

(23) 
(24) 

(25) 
(26) 
(27) 

(28) 
(29) 

(30) 
(31) 

(32) 

(33) 
(34) 



Equations JBJ— ©, (fT9)l . ((2T1) . ([22)1. (I3T1) . ((32]) and (|34)l yield exclusive expression for pressure (p), cosmological 
constant (A), Gravitational constant (G) and bulk viscous stress (II) as follows, 



P = 



3j(2n + l)fci 



87rG ^3m(n + 2) 2 (fci/j + k )^ 



A 



(2n + l)fci 



(n + 2) 5 



3^— 



(fcif + fco) 2 m(k 1 t + k ) 



G = G [£(kit + k }^ 



n 



(2n+ l)fci 



87rGo(n + 2) 2 [£(fcit + fc )]^ 



2fci 



3(1-7)£ I2 5- 

m(k 1 t + k ) (fci/j + fco) 2 



(35) 

(36) 
(37) 
(38) 
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Figure 1: Plot of Gravitational constant (G) versus cosmic time (f). 



We observe that model has singularity at t = — % which can be shifted to t = 0, by choosing ko = 0. This 
singularity is of point type as all scale factors vanish at t — — %. The parameter p, p and A start off with 
extremely large values. From (|28[). it can be seen that the spatial volume is zero at t = — % and it increases with 
cosmic time. The parameter H x , H y , H Zl H, 9 and a 2 diverse at initial singularity. These parameters decrease 
with evolution of Universe and finally drop to zero at late time. Fig. 1 depicts the variation of gravitational 
constant G versus time. From we observe that A(t) is decreasing function of time and A > for all times. 
Fig. 2 shows this behaviour of cosmological constant A(t). Thus the nature of A in our derived model of the 
Universe is consistent with recent SN la observations. 



In EIT, the bulk viscous stress II satisfies a transport equation given by 



n + rIL 



-3£ff 



-Til 



3H 



T 



T 



(39) 



where, r is the relaxation coefficient of the transient bulk viscous effects and T > is the absolute temperature 
of the Universe. The parameter e takes the value or 1. Here e = 0, represents truncated Israel- Stewart theory 
and e = l, represents full Isreal-Stewart (FIS) causal theory. One recovers the non-causal Eckart theory for 
t = 0. 

Maartens (1995) has pointed out that the Gibb's integrability condition suggest if the equation of state for 
pressure is barotropic (i. e. p = p{p)) then the equation of state for temperature should be barotropic (i. e. 
T = T(p)) and it may be expressed as 

dp(p) 



Toe 



From equations (131)) and (|34j) . we obtain 



p + p{p) 
T = T oP ^ , 



(40) 
(41) 



where Tq stands for a constant. 

Using Eg. (1341) into Eq. (|41[) . we obtain the expression for temperature (T) in terms of cosmic time (t) as 



T = Tn 



3(2n + l)fe! 



8irGo£zm(n + 2) 2 (M + ko) 1 



(42) 



From equation (1421) . it is evident that temperature is decreasing function of time. The variation of temperature 
versus cosmic time for 7=5 (radiation dominated era) and 7 = 1 (stiff fluid dominated era) has been graphed 
in Fig. 3. It is clear that temperature of Universe decreases sharply for stiff fluid and approaches to small 
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Figure 2: Plot of Cosmological constant (A) versus cosmic time (t). 




G 



positive value at late time, as expected. 

Bulk Viscosity in Eckart's Theory: The evolution equation (f3"9")l for bulk viscosity in non-causal Eckart's 
theory reduces to 

n = -3£H . (43) 

With help of equations ([25jl , (|38|) and (|43|) , we have the relation between bulk viscosity coefficient (£ ) and cosmic 
time (t) as 



(2n + l) 



87rG (n + 2) 2 [^(A:it + A:o)] T 



3(1 - j)£— 



2k x 



(fat + h ) 



(44) 



Bulk Viscosity in Truncated Theory: ft has been already pointed out that in truncated theory (i. e. e = 0), 
the evolution equation (13"9")l for bulk viscosity reduces to 

II + rn 4 = -3£H . (45) 

Following, Singh et al (2009), the relation between r and coefficient of bulk viscosity £ is given by 



(46) 



This relation is physically viable because the viscosity signals do not exceed the speed of light. Thus the equation 
(1351 leads to 

(47) 



n + -n 4 = -3£H . 
p 

Using equations (f25j). (l34l) and (|38|) into equation (|47l) . we obtain 



3(1 - tK 1 ^- - 2fe 1 (fe 1 f + fc )m 



[e(fat + k )f][fak 3 (k 1 t + k ) - k 4 } ' 



(48) 



where 

. _ (2n+l)fci 



87rG m(n+2) :2 ' 

_ ( m+ 3)(l- 7 )+lt 
^3 — j 



u 1 



k A - — 



Bulk Viscosity in FIS Causal Theory: Using equations (1531) and (13"9"|) . the transport equation (l3"9l reduces 
to 



n + in 4 = -3ff£- |- 

p 2p 



3H 



(1 + 2 7 )P4 



(49) 



(1+7)P 

Further, Using equations (|25l) . (IM)) and (|3"51) into equation (|49l) . one can easily obtain the relation between bulk 
viscosity coefficient (£) and cosmic time (t) as 



mfakzildt + k 2 ) 2 ^(t) 



[t(fat + k )]ir [fak 3 (fat + k ) + k 5 (fat + k ) 2 ip{t) - fc 4 



(50) 



where 



m[3(l+7)-(l+2 7 )(m+3)] 



3(1-7)^^ 

m(kit+k ) (fcit+fco) 2 



2fci 



4 Distance Modulus Curves 

The distance modulus is given by 

H = 5 logd L + 25, (51) 
where is the luminosity distance and it is defined as 

d L =r 1 (l + z)a , (52) 
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where z and ao represent red shift parameter and present scale factor respectively. 

For determination of r\, we assume that a photon emitted by a source with co-ordinate r — r\ and t = t x and 
received at a time tj by an observer located at r = 0. Then we determine t\ from 



dt 

n = I - 



to , 

I, " 



(53) 



Equation (j3"Tj) can be rewritten as 

a = l^i , (54) 

where t = k±t + ko. 

Solving equations (|5 1 [) — (|54[) . one can easily obtain the expression for distance modulus (p) in term of red shift 
parameter (z) as 

fa , n2 



/i = 5 log 



2H (1 + z, 



((1 + -) 2 



1 



25 



(55) 



Table: 1 



Redshift (z) 


Supernovae la {^l) 


Our model (p) 


0.014 


33.73 


33.81 


0.026 


35.62 


35.17 


0.036 


36.39 


35.89 


0.040 


36.38 


36.13 


0.050 


37.08 


36.63 


0.063 


37.67 


37.14 


0.079 


37.94 


37.66 


0.088 


38.07 


37.90 


0.101 


38.73 


38.22 


0.160 


39.08 


39.29 


0.240 


40.68 


40.26 


0.380 


42.02 


41.40 


0.430 


42.33 


41.71 


0.480 


42.37 


42.01 


0.620 


43.11 


42.67 


0.740 


43.35 


43.15 


0.778 


43.81 


43.28 


0.828 


43.59 


43.46 


0.886 


43.91 


43.64 


0.910 


44.44 


43.72 


0.930 


44.61 


43.78 


0.949 


43.99 


43.83 


0.970 


44.13 


43.89 


0.983 


44.13 


43.93 


1.056 


44.25 


44.13 


1.190 


44.19 


44.47 


1.305 


44.51 


44.73 


1.340 


44.92 


44.81 


1.551 


45.07 


45.235 



The observed value of distance modulus (ijl(z)) at different redshift parameters (z) given in table 1 (SN la 
Data) are employed to draw the curve corresponding to the calculate value of n(z). Fig. 4 shows the plot of 
observed /i(z) (dotted line) and calculated /i(z) (solid line) versus redshift parameters (z). 

5 Concluding Remarks 

In this paper, we have presented exact solution of Einstein's field equations with variable G and A in LRS 
Bianchi type I space-time in presence of imperfect fluid. The main features of the work are as follows: 



8 



46 r 



44 - 



40 ■ 



38 ■ 



36 ■ 



34 - 



32 



0.2 0.4 0.6 



1.2 1.4 



1.6 



Figure 4: Plot of distance modulus (/i) versus redshift (z) for Supernova data (dotted line) and for our model 
(solid line). 



• The derived model represents the power law solution which is different from other author's solution, 
seems to describe the dynamics of Universe from big bang to present epoch. 



It 



• The cosmological constant (A) is found to be decreasing function of time and it approaches to small pos- 
itive value at late time. A positive value of A corresponds to negative effective mass density (repulsion). 
Hence we expect that in the Universe with the positive value of A, the expansion will tends to accelerate. 
Thus the derived model predicts accelerating Universe at present epoch. This is in the favour of recent 
supernovae la observations. 

• The temperature of Universe in derived model is infinitely high at early stage of evolution of Universe 
but it approaches to small positive value at later stage. This means that temperature is also decreasing 
function of time. The same is predicted by CMBR observations. 

• If we choose n = 1, the mean anisotropy parameter vanishes. Therefore isotropy is achieved in the de- 
rived model for n = 1. Also we see that for £ = 1 and n = 1, the directional scale factors vary as 
A(t) = B(t) = a(t), therefore metric ([1]) reduces to the fiat FRW space-time. Thus I — 1 and n = 1, turn 
out to be the condition of flatness in the derived model. It is important to note here that for n = 1, shear 
scalar vanishes but the bulk viscosity contributes to the expansion of Universe and for positive value of n, 
the bulk viscosity coefficient (£) decreases with time. 

• The distance modulus curve of derived model is in good agreement with SN la data (see Fig. 4 and Table 
!)• 



The age of Universe is given by 



T - 1 H- 1 k ° 



Finally, the model presented in this paper is accelerating, shearing and starts expanding with big bang singu- 
larity. This singularity is of point type singularity. 
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